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Abstract 

We develop a 5d model for Electro Weak physics based on a non compact warped extra 
dimension of finite length, known as the soft wall scenario, where all the dynamical degrees 
of freedom propagate in the 5d bulk. We solve the equations of motion and find the allowed 
spectra, showing that the mass of the lightest fermionic mode behaves as a power law of 
the effective 4d Yukawa coupling constant, with the exponent being the corresponding 
fermionic 5d bulk mass. Precisely this non universal behavior allows us to reproduce the 
hierarchy between the Standard Model (SM) fermion masses (from neutrinos to the top 
quark) with non-hierarchycal fermionic bulk masses. 
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1 Introduction 



The existence of warped extra dimensions has been a source of prolific investigations for 
the last 10 years because it was shown that it could provide with an alternative explanation 
of the hierarchy between the Planck and the Electro Weak scales [1, 2]. The pioneer work 
by Randall and Sundrum [1] considers a slice of a five dimensional anti de Sitter space 
spanning between two 3-branes, the UV one accounts for the gravity (Planckian) scale 
and the IR brane is chosen to have a scale of the order of the ElectroWeak breaking by 
the exponential suppression of the AdS metric. Another aspect of this set-up that has 
attracted a lot of attention is its connection to the AdS/CFT duality [3] which states 
that this 5d theory is equivalent to a 4d CFT coupled to gravity and the interpretation 
of the IR scale is that it corresponds to the breaking of that conformal symmetry and the 
appearance of a mass gap and resonances [4] . 

More recently, an alternative scenario with a (non compact) warped extra dimension 
has been proposed [5, 6, 7], the so-called soft wall as opposed to the hard wall termination 
of the Randall-Sundrum model. The motivation for the new scenario, as pointed out in 
reference [5] , comes from the fact that a suitable soft wall configuration may lead to linear 
Kaluza- Klein (KK) excitations (m n ~ n), unlike the quadratic excitations found in the 
hard wall models, and hence it may provide a description of the linear spectrum of mesons 
in QCD [8]. This is what has been called the AdS/QCD approach. 

In the soft wall case the (effective) metric is not AdS, instead the warped factor decays 
faster enough to make the extra dimension of finite length and in this sense a virtual 
IR boundary is located at infinity but with a finite distance from the UV one and the 
departure from the AdS behavior is associated to a smooth VEV acquired by some dilaton 
field, however, since no real IR boundary exists all the matter/radiation content propagate 
in the 5d bulk, thus the suppression due to the metric is universal (and is absorbed through 
the wave function normalization). The possible suppression on the lowest modes comes 
from the binding of the Schrodinger-like potentials governing the profiles of the wave 
functions [6, 7]. 

In this work we try to reproduce the hierarchy of the fermion masses of the Standard 
Model embedding it on a soft wall background. It is worth remarking, nevertheless, that 
our main concern here is to reproduce the correct order of magnitude for the masses, a 
more realistic study taking into account all the families and all the ElectroWeak constrains 
would need to incorporate the CKM matrices, which we leave for future investigations. In 
any case, the incorporation of the mixing angles would not change the main results we are 
presenting. Finally, since the topic is relatively new and relatively absent in the literature 
we think a detailed derivation, although overlapping previous publications, is worthwhile. 

The paper is structured as follows: In section 2 we present the model together with the 
notations and conventions. Section 3 is devoted to solve the equations of motion for the 
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fermion, Higgs and gauge sectors. Section 4 applies the previous results to the hierarchy of 
the fermion masses and section 5 contains the conclusions. Finally, at the end of the paper, 
we have included appendices developing some technical aspects of particular calculations. 




2 The Model 

Before presenting the model we will set up the notation used: Capital letters from the 
beginning of the alphabet (A, B, ■ ■ ■ ) represent (local) Lorentz 5d indices while those from 
the middle and end of the alphabet (M, N, R,---) stand for general coordinate 5d indices. 
Finally, greek letters take into account 4d Lorentz indices. 

The 5d Lorentz representation is spanned by the gamma matrices j A = (7^,— i7 5 ), 
with 7^ the usual 4d gamma matrices and 

7 5 = 

Our model is defined in a five dimensional manifold with (UV) boundary, £ = R 4 x X, 
with X = [zq, 00) and z$ > 0, assuming also zq <C (Mvk)" 1 being My/ the Electro Weak 
scale. The background metric is AdS, that is 

9 = f(z) V5 , (1) 

where f{z) = k/z, k being the radius of curvature and as a matter of fact we will take 
zq ~ k. In addition 775 stands for the 5d Minkowski metric with mostly negative signature, 
namely 775 ab = diagonal (+1, —1, —1, —1, —1) and z represents the fifth coordinate. The 
propagating degrees of freedom consist of an SU(2) l doublet of 5d Dirac fermions: and 
two singlets an SU(2)l doublet of complex scalars, H\ and a gauge field, Am, lying 
in the adjoint representation 1 of SU(2)l x U(1)y, where the lightest modes of the above 
fields are to be identified with the SM fermions, Higgs and Electro Weak gauge filelds, 
respectively. For the 5d fields, the bulk-brane action is 2 



-n 7 M v L M ^ L + -n 7 M ^ n + m + m 



+ (Ai eijVi&R H*i + A 2 H**^ + h.c.) 

+g MN (D M H)1 D N H - ml^H + -L g MR g NS tr {FmnFrs} 
Jg-e-* \\ Q k 2 (\H\ 2 - v 2 ) 2 - el 1 - - 



(2) 



x To avoid important quantum corrections to the p parameter the gauge group should be SU(2)l x 
SU (2)rxU(1)y [6, 7] although for the sake of simplicity we take the usual gauge group. For the calculations 
of the spectra this technicality does not affect. 

2 The presence of the (5,5) component of the inverse vielbein, e|, in the brane action is not a problem 
because the (effective) background metric itself breaks the isometry group to 4d Lorentz-Poincare. 
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where Fmn = d[ M A N ] + L4m,Aw] and Am = g$ A S M T S T s being the group generators 
verifying [T s , Tt] = i/Jt?r and 55 the (5d) gauge coupling constant, A 1,2 are the 5d Yukawa 
coupling constants, M is a (constant) bulk mass and m 2 = ^ —2a^ z 2 with a, Ao being 
dimensionless constants whereas vq has dimension of (energy) 3 ^ 2 . v L j is the Levi-Civita 
tensor and 4> is the dilaton field. In addition 

D M = d M + iA M , (3) 
V L M = D M + E < M T£ B r lAD X BD , (4) 
V M = d M + E M T£ BVAD Z BD , (5) 

are gauge and Lorentz covariant derivatives, where E M is the vielbein field, verifying 
9m N = E M E B rjsAB- The gamma matrices depend on the spacetime coordinates according 
to 7 M = e^j A with being the inverse of the vielbein field verifying qm n e A e B = Vsab, 
i.e. they are the (local coordinate) components of orthonormal vector fields (e^ = c^Om)- 
Furthermore, T, AB = | [7 j4 ,7 B ] are the Lorentz group generators in the representation 

diagonal (SL(2,C) x 5L(2,C) t ^ 1 ) , 

and the symbols verify V es ec* = with V the Levi-Civita connection compati- 

ble with g. The relation with the Christoffel symbols, T MN , is then T^ c = e B EfjdM&c + 
e M e^,E^T MN . In the case of the metric (1) the Lorentz covariant derivative is 

^ = d,+ l -l"ljj, (6) 

V? = d 5 . (7) 

For the dilaton field we will assume a z-dependent VEV. The particular bulk-brane 
gravitational dynamic leading to that profile for the dilaton is beyond the scope of the 
present work, although in Ref. [7] an underlying gravity model leading to a soft wall 
configuration can be found. Here we will take ((/)} = /j, 2 z 2 with \i some mass scale. As for 
the Higgs bulk mass is concerned, its z dependence can be understood as a consequence 
of the coupling to the gravity sector. We will also not address the question of the different 
scales in the problem but supposed that their stability is granted by an underlaying UV 
theory. 

Finally, we have taken the same (z-independent) bulk mass for Left- and Right-handed 
fermions for the sake of simplicity at solving the equations of motion. Different masses, as 
it can be seen in the next section, substantially increase the difficulty as it is not possible 
to decouple easily the first order fermionic equations of motion. 

3 Solving the equations of motion 

We will start by finding the VEV of the Higgs following the derivation found in Ref. [7] . 
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3.1 Higgs VEV 

Consider the purely scalar sector 3 out of the action (2) and assume for the moment 
H(x, z) = ft(z), then, the minimal action principle yields the equations of motion 

j -</) ft'-f 2 m 2 h ft = 0, (8) 

together with the boundary conditions 

ft' -2X k 2 (ft 2 -v 2 )ft\ zo = 0. (9) 

Redefining the scalar as ft = z a H we find the equation 

w ,/ + ^2a-3_ 2 ^ w , = ()j (1Q) 

and thus the general solution is 

ft = z a (c + c j* dc C 3 - 2a e" 2 ^ , (11) 

although the normalizable solution is obtained with c = and hence the boundary condi- 
tions set the constraint 

c[a-2\ Q k 2 {c 2 zl a -vl)zv] = 0, (12) 

whose solution is 

' c = 

(13) 

2 _ I 1 1 / a 1 2 



C - [zoj \2z \ k* +V < 

Now if we compute the energy density per unit 4d volume [9] for the solutions to the 
system (8)- (9) we find 

Ssi = I v^e-^ (-g 55 d 5 ftd 5 ft + m\ft 2 ) + / ^e-*\ k 2 (ft 2 - v 2 ) 2 

Jt JdT 

= J^V9 e-^ft [g 55 V 5 d 5 ft - <? 55 d 5 8 5 ft + m 2 h ft] 
+ f \ftft' + X k 2 (ft 2 - vlf 

JdT 1 i 

=X k 2 [ ^ge-*(v 4 -ft 4 ) . (14) 

JdT 

The energy density per unit 4d volume is then decreasing around the trivial solution, 
ft = 0, for Ao > and therefore in this case the non trivial solution is the absolute 
minimum. 

Now that the Higgs profile has been set, we then proceed by solving the equations of 
motion for the fermion sector according to three polynomial cases for the Higgs VEV. 



'Wo consider the Higgs doublet in the unitary Range. Mini is // = [ ^ J with '0 
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3.2 Fermionic sector 

Applying the variational principle to the fermionic sector 4 of the action (2) one finds a 
bulk piece whose vanishing induces the following equations of motion 

ie^ 7 A ^n " 2 47 + Mn + H + ^ K = , (15) 

ie^^V^ - \e 5 Al A 4>^R + M* 1 R + Aie^tf = , (16) 

ief^V^l - \eh A ^ 2 R + M* R + \ 2 tf n = (17) 
while the boundary variation yields the constraint 

el, (l - <x 3 ® 7 5 ) ( 11 ) =0. (18) 
V R J z 

where 03 is the z-Pauli matrix acting on (L, R) space. If we refer and ^ R as simply 
and ^r, respectively, the equations of motion and the boundary conditions adopt the 
generic form 

ief^V^L - i4 7 V*L + M * L + \S}(z) * R , (19) 
ief^V^R - X -e\ 7V** + MV R + \?>(z)y L = 0, (20) 

e|(l-a 3 ®7 5 ) ( l L ) =0. (21) 

\ / ZQ 



and upon the redefinitions 5 



*L,R=^r 2 e^(iP + Tip-) , (22) 



we find 

\^ ±+1 b d^ ± + {M±m D )ftlj ± = 0, (23) 
with mo = Xfy^z), while the boundary conditions take the form 

(l-^7 5 )(^ + )=0. (24) 
Considering now an orbifold-like decomposition 6 , that is 

** - ( h if ±{ <^ ) ■ < 25 > 

V 9±{z)x±{x> x ) / 



4 Considering the Lorentz covariant derivative only. 

5 With different L and R fermionic bulk masses one can not decouple the system with a global rotation. 
6 Notice that the normalization of the kinetic term in these variables is simply \h\ 2 . 
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where £±, x± are Weyl spinors verifying \<J^d^± = m\±, m being the physical eigenmass, 
the first order coupled differential equations (23) can be written as 



h'± + (M± m D ) fh± = -mg± , 
g'±-(M± m D ) fg± = mh±, 



(26) 
(27) 



and differentiating once more they transform into two decoupled 2nd order differential 
equations 

til 



± 



n 
9± 



f (M ± m D f -f'(M± m D ) T f m' D - m 2 \ h± = , (28) 
f (M ± m D f + f(M± m D ) ±fm' D - m 2 ] g± = . (29) 

In addition the boundary conditions (24), in terms of the new variables, translate into 7 



(1-ai) 



(l + ^i) 



5+ 
9- 



0. 



(30) 



In the following we explicitly solve these equations for three polynomial behaviors of the 
Higgs VEV. 

• Constant profile (a = 0): Sj(z) = vq 

In this case the equations take the form 



d a± 



(31) 
(32) 



de 2 

d^ 2 



9± 



e 2 

a± (g± - 1] 
9 2 



g± + g± = o, 



(33) 
(34) 



with a± = Mk ± Xvok and 9 = mz. Redefining the fuctions as h± = Vdi)± and g± = 
\[Q g ± we find the Bessel equations 



d 

d# yaH 

d / d 



d6 V d6 



t>±) + 
0± ) + 



a± + 



- \OL± ~ 



f)± = 0, 



0± = O, 



(35) 
(36) 



7 Eq. (24) implies \+ = K X- with k a constant (and an analogous relation for £±) however by redefining 
the Weyl spinors one can absorb the constants. 
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and, accordingly, the solutions are 

rJ 



h± = Ve 
g± = Ve 



Hi .1 [ a± + \,o) +HlY U± + \ 



G J ± J[a±- 



1 



+ G±Y [a±- 



1 

2'' 



for a± ± i non-integer. One can easily check that 



d0 



+ a± 



V~9Y(l 



and hence the first order constraints (31)-(32) impose 



G 



± 



-H. 



± > 



Gl = -Hl. 



(37) 
(38) 

(39) 
(40) 

(41) 
(42) 



Since we have enough degrees of freedom to solve the boundary conditions (30), no further 
restriction on the mass eigenvalue is found. Thus, the spectrum is continuum and there 
is no mass gap. 

• Linear profile (a = 1): fi(z) 

The equations of motion are now 



\k 3 . 



de 



/ Mk A ft 



TTi 



9± 



Mk A fl\ m 
^ ± W = kl) g± = ^ h± 



(43) 
(44) 



de 2 
d 2 



h± 



>9± 



Mk {Mk + 1) , MXp,y/k 1 

e 2 7 e + 1 

Mk{Mk-l) MXfiVk 1 



Q2 



7 e 



h± = 0, 

g± = o, 



de 2 

with = 2jz and 7 = \J^rfi 2 — m 2 . Then with the redefinitions 

h± = e- / 2 e 1+Mfc f,± , g± = e- @ / 2 Q 1 -^ B± , 



(45) 
(46) 

(47) 



Eqs. (45)-(46) turn into the following confluent hypergeometric equations (alsol known as 
Kummer equations) 



^2 &± + ( 2 + 2a ~ ) ^ &± " (! + « ± $ &± = , 
e^ 0± + (2-2a-e)^ B± -(l-a±/3) 0± =0, 



(48) 
(49) 
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where a = Mk and = MA ^ V ^ . For 2 ± 2a ^ Z the general solution is (see appendix A) 

f)± = #± 1F1 (1 + a ± 0, 2 + 2a, 9) + H± U (1 + a ± 0, 2 + 2a, 9) , (50) 
0± = G± iFi (1 - a ± 0, 2 - 2a, 9) + G± £/ (1 - a ± /3, 2 - 2a, 9) . (51) 

Furthermore, one easily show that 
d 



J |9 1+a e- / 2 1 F 1 (l + a±/3,2 + 2a,9)} 
(1 + 2a) 9 a e~ e/2 1F1 (-a ± 0, -2a, 9) , 



d9 V 9 2a 



(52) 



_d_ 
d9 



+ i . ^ |9 1+a e~ 0/2 C/(l + a±/3,2 + 2a,9)} = 



" ± JL 

9 2a 



Q a e -e/2 v ( _ a ± /?) _ 2aj e) 



(53) 



and so for M — > — M and /i — > 5, thus the general solution to the coupled system (43)-(44) 
is given by 



h± = Q 1+a e~ 0/2 [/± 1F1 (1 + a ± 0, 2 + 2a, 9) 
+ u±U (1 + a± 0,2 + 2a,@)] , 



(54) 



7 (a T /?) 



,/± 27 (1 + 2a) 1^1 (-a ± ft -2a, 9) 



m 



+ u± 



m a 



C/(-a±/3,-2a,9) 



As far as the spectrum is concerned we find two distinct regions: 



1. m 2 > 



A 2 /t 2 
k ' 



i.e. 9 = i |9| (continuum spectrum) 



(55) 



In this case none of the above linearly independent solutions are normalizable and 
analogously to the previous case the boundary conditions at zq are satisfied without 
further restrictions on the physical mass. 

\ 2 ~ 2 

2. m 2 < —jjf- (discrete spectrum) 

For real 9 only the [/-type solution is normalizable since i-Fl (a,0, 9) ~ Q a ~P e® 
while U (a, 0, 9) ~ Q~ a as 9 — > 00. We then set f± = and thus a non trivial 
solution to the boundary conditions (30) requires the vanishing of the function 

(a + 0) U (1 + a + 0, 2 + 2a, 9 ) 17 (-a - /3, -2a, 9 ) 
+ {a-0) ?7(l + a-/3,2 + 2a,9o)?7(-a + /3,-2a,9o) , 
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with 9 = 272:0 < 1. 

Therefore the linear case predicts a finite set of eigenmasses followed by a continuum of 
eigenvalues. However the mass restriction does not allow a light mode. 

• Quadratic profile (a = 2): S)(z) = \ j£ with p, = (^Y J^j-+v%k 3 . 



In this last case, the (first order) equations of motion turn into 

-^h± + ^(a±()h± = -5-^=g±, (56) 

^9±-^(<x±Qg± = 8-j=h±, (57) 

with C = (j^j Co, a = Mk, 5 = and Co = + ""o^ 3 - Redefining the functions 



as 



h + = c *e 2 c f)+, 5- = C 2 e 5C B _, 

1 + Q 1 1 — Q 1 > 

fc- = C 2 e 2 c t)_, 5+ = C 2 e 2 C 0+ , 



the second order differential equations become the following Rummer's equations 

C^+Q-«-c)f/+ + ^f) + = 0, (58) 

C fl ':+Q + a-c)fl'_ + jfl- = 0, (59) 

C f)- + (I + a " C) f>- - (l - j) f)- = , (60) 



C0'; + g-«-c) 9 ' + -(i- j) 0+ = o, 

such that the (normalizable) solutions are 

Hrrf ^ 1 A 

h+ = a+C *e *tU\--,--a,t) , 

l+a 1, / (5 2 3 \ 

/»- = «- C~ e~2 C ?7 f 1 - -, - + a, CJ , 



(61) 

(62) 
(63) 



a,C , (64) 



5 _ = fe_ C f 17 f-^, 1 + a, cV (65) 



while the first order constraints impose 



b + = -^5a + , a_ = ^5b_, (66) 
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Figure 1: Normalized fermionic mass quantization condition as a function of S for Co = 
10 -3 and a = 3 (solid line) and a = 5 (dashed line). 

then the solutions are given by 

I r*e-' c tf (-?,$- «,<)<(*) \ .... 
* += Uc^eV(i-?,t- a ,c)x(*)j' <67) 

/fC ! Pe-lCv( l -?,i + a,c)«.)\ 
*- = i <Se-l^(-f,i + a ,c)x W J" (68) 

Finally, a non trivial solution out of the boundary conditions (30) set the restriction 
/ 5 2 1 \ / 5 2 1 \ 

_^2 Co[/ A_^3_^ \ ^^^^^ \ =Q 

This condition quantizes the possible physical masses and indeed it presents a mass gap. 
To see this notice that the confluent hypergeometric function U (a, b, z), as a function of a, 
oscillates for negative values of this parameter, then we expect to find oscillatory behavior 
for 5 > 2. Actually, in Fig. 1 we plot the above mass quantization condition for £o ~ 10 -3 
and two different values of a. 

Further assuming that Co "C 1 (at the end of section 4 we justify this smallness) and 
for \a\ > 1/2 we find a light smallest Dirac eigenmass, given by 

m 2 k 2 *2 ] (Co) H+ ^ • (69) 

For \a\ < 1/2 the lowest eigenvalue is not power law suppressed [7]. Indeed the authors 
of Ref. [7] propose such different behavior as a source of the hierarchy. Nevertheless, as 
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we shall see, the hierarchy can be induced within the power law behavior. As for the 
spectrum is concerned, we are interested in a discrete set of KK excitations, hence from 
now on we will consider a quadratic profile for the Higgs VEV and as pointed out in the 
previous section, the lightest modes will be identified as the SM fields. 
To finish this section we will present the solutions for the bosonic sector. 

3.3 Bosonic sector 

3.3.1 Electro Weak equations of motion 

The variational principle applied on the gauge piece of the action (2) (taking into account 
the Higgs VEV) sets the equations of motion 

±f N {V9 9 MR 9 NS d [R A s] } - 9 MN 9 RS d [N A R] d s <t> - 9^ 2 g MN A N = , (70) 

along with the boundary conditions 

F 5 »=g 5M g» N d [M A N] \ zo =0. (71) 

Applied to the AdS metric (1) in the gauge ^.5 = and 8m A M = we find the equations 
of motion 



(-□ + d 2 5 + -(f/ds- <?|^ 2 / 2 ) ^ = , 



(72) 



d 5 A\ Q = 0. (73) 
Decomposing the gauge field as A^(x,z) = v lx {x)g{z) with Qf^ 1 = — m 2 v M we find 



g" + p - <t>'g' - (glrff - m 2 ) g = . (74) 



Then with the redefinitions 



g = f-W e** e-* e 9~* g 

a _ /-1 1 g'jp 2 „2„2 



(75) 



we find 



'^•-»s» +m ' , = - <76) 



where m = ^ and 7 = 1 + . Finding the general solution to this particular case of 
Kummer differential equation is somehow involved. Here we simply give the normalizable 
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solution, although a detailed derivation can be found in appendix A. For the gauge field 
the rn-th mode is 

A"{x, z) = v»(x) e^ 1 '^ 2 * 2 g ( 7/ uV) , 

with 



where 



g (9) = 1 - m 2 9 ln6 xFx (l - m 2 , 2, 9) + f ^^y u (*) > 



00 r fi 2 + 

u ( g ) = E Lr/^M [^o(l + fc) + ^o(2 + fc)-^ (l-m 2 + fc)](9* 



fc=0 



A:!r(2 + fc) 



and ( x ) = is the digamma function. The boundary conditions (73) are 

^(l-7)0(#o)+70'(#o) = O. 
Assuming then that 9$ "C 1 we find for the lowest mode (m 2 <C 1) 

1 - 7 - 2 7 m 2 [1 - (2) + hi0 o ] + O (9 ln9 ) = , 

and hence 

1 

m 2 w » 2^7(7 - 1) 



Mol 



(77) 
(78) 



(79) 



(80) 



(81) 



Notice that 7 is always greater or equal than 1, and thus the mass is always real. On 
the other hand, when the coupling to the Higgs VEV is switched off (55 — > 0), which is 
the case of the photon, the lowest eigenmass is zero. Indeed, for 7 = 1 and m = the 
equations of motion can be written as 



fe-<t>d 5 (fe-tdaA") =0, 



(82) 



and together with the boundary condition fe ^ d^,A^\ Z0 = 0, both impose a constant 
profile for the massless mode. 



3.3.2 Higgs equations of motion 

As far as the Higgs is concerned, here we will analyze the perturbations around the back- 
ground found in section 3.1, that is H(x, z) = Sj(z)+H(x, z), which gives the free equations 
of motion 

(83) 



UH - d 2 H - (3^ - ^ 8 5 H + / 2 m 2 H = 0. 



with OH = —m 2 H, while the boundary variation imposes the vanishing of the expression 



d 5 H - 2X k 2 



9) + H 



(84) 
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at the boundary z = zq. For solving the equations of motion we redefine the the function 
as H = / _3//2 r4f) with r = /j, 2 z 2 and then we find the Kummer equation 

Again we refer to appendix A to find the details of the solution and here we simply write 
down the normalizable one, which turns out to be 

2 \ / 2 

m \ , „ / m 



H(x,r) = H(x)T[ -— 2 )r\nr 1 F 1 -^.l.r 



4fi 2 J \ 4/j, 2 

"V, (86) 



- 4r(W r (l + fc) - *° (- v + fc ) 



or shortly: H(x, z) = H(x) 7] (/x 2 z 2 ). Assuming that To <C 1 the boundary conditions can 
be expressed as 

1 + lnr - 2* (1) + ) \( + O (r lnr ) = , (87) 
and since ~ — \ for e <C 1 we find that the lowest Higgs mass is 

"* as4 "W < 88 > 

4 Phenomenological predictions on the hierarchy of the SM 
fermion masses 

In this section we will apply the previous results to the fermion masses. In particular we 
will use them as a possible source for the hierarchy of the masses. 

Whenever we have normalizable solutions, a 4d effective action is well defined. Starting 
from the action (2), by means of partial integration and upon the redefinitions (22) we 
find the action (in appendix C we give a proof for the orthogonality of the normalizable 
solutions) 



n>0 

where 8 



J d 4 ^ (iVfW f (n) - mj n) vfV f (n) ) , (89) 



/ /(") 

^ -Us 



Recall that the boundary conditions (24) set £ + = £_ and \+ = X- an d thus we are left with a single 
Dirac fermion. 
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with f labeling the fermion flavour. In addition we will assume zq ~ k. 

Notice, however, that the values of the 5d couplings are physically meaningless by 
themselves, instead one should consider the effective 4d couplings obtained upon the inte- 
gration over the fifth dimension, and for this aim we need to normalize the wave functions. 
The normalization constants are given by the expressions 



N 
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for the fermions, while for the Higgs one finds 
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As for the effective Yukawa coupling is concerned, we can read it off from the action (2) 
^e-'t'XH^L^R + ^R^L) = \fH$+1>+-$-il>-) , (93) 
and therefore we find the coupling 

Y (n,m,r) = y M R{n) {%) ^M (x) ^M (x) + ^ > (g4) 

with 
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(95) 



where, for simplicity, we have omitted the flavor index in the integrals. Thus the effective 
SM Yukawa coupling would be the one between the lightest modes, which will be referred 
to as simply y^ d = y^'°^ ■ In appendix B we give an estimation of the previous hyper- 
geometric integrals based on the smallness of the parameters tq , Co and the exponential 
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-log |0 (mk) 




Figure 2: — log 10 (mk) as a function of a for k 1 ~ W 19 GeV and either y^ < 1 (dashed 
line) or y^ ~ 0.1 (solid line). 

suppression of the integrands. According to that approximation we can estimate the order 
of magnitude of the fermionic physical (Dirac) mass as 

< 

Recall that at = Mtk, Mt being the 5d fermion bulk mass. 

We will be interested, nevertheless, in the case \a\ > |. Indeed this power law de- 
pendence on the Yukawa coupling allows us to reproduce the hierarchy between the SM 
fermion masses with a tiny variation in the 5d bulk masses. 

To extract some numbers we will take the scale mass ji associated with the dilaton 
VEV as fjL ~TeV, indeed this is needed for the gauge and the Higgs masses, Eqs. (81) 
and (88), respectively, to be at the Electro Weak scale, and similarly we will take the 5d 
VEV of the Higgs as vq ~ ^TeV and |Ao| < 1. For the curvature, however, we will 
consider two extreme situations: 

1. Planckian curvature, k^ 1 ~ 10 19 GeV. 

As we can see in Fig. 2 with order 1 effective Yukawa couplings, y^ < 1, we find the 
correct order of magnitude for the top quark mass (mtk ~ 10~ 17 ) for \at\ — 13 and 
for the neutrinos masses (m u k ~ 10~ 28 , 1CP 31 -depending on wether they are light, 
~ eV, or ultralight ~ meV-) with \a u \ — 20, 22, respectively. 

Of course, the smallness of the mass scale fi compared to the AdS curvature, fj,k ~ 
10 -16 , and the stabilization of the dilaton VEV deserve an explanation that we are 
not addressing in this work. 
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Figure 3: — log 10 (mk) as a function of a for k 1 ~ W^GeV and either y^d % 1 (dashed 
line) or y&d ~ 0.1 (solid line). 

If we, instead, consider slightly suppressed 4d Yukawa couplings, y^d ~ 0.1, the 
suitable values of a are: \at\ — 10 and \a u \ ~ 15, 17. 

2. LHC curvature, Ar 1 ~ 10 4 GeV. 

In this case (as Fig. 3 shows) with 0(1) Yukawa couplings we reach the order of 
magnitude for the top quark (mtk ~ 10~ 2 ) with \at\ — 5 while for the neutrinos 
(m u k ~ 10 -13 , 10~ 16 ) we have to consider \a u \ — 18,21, respectively. 

And as before, taking y^d ~ 0.1 we instead find \at\ — 2 and \a u \ ~ 12, 14. 
Concerning the parameters Co) T 0) from Eq. (69) we can express Co as 



One should compare the previous results with the case of Randall-Sundrum, where 
the hierarchy of fermions, at least if one forgets about neutrinos, is obtained for 0(1) 
parameters whereas here we need to have 0(10) parameters. Although this is not a big 
difference, the reason why we have such a result is that we are enforcing both the left- 
handed doublet and right-handed singlet to have the same 5d mass whereas in the case 




from where we find 



< 
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of RS the prescription taken is cl = 1/2 — cr. Have we taken a different prescription we 
would have obtained different results, but we would have lost any possibility of solving 
the problem analytically. Since the whole point of this paper is just to show how one can 
generate the hierarchy and not to give a complete model we have preferred to be as simple 
as possible. 

Let us finished this section by saying that all the masses that we generate in this 
set-up are Dirac since they come from a Yukawa coupling of the Higgs with two different 
fermions. This will imply that neutrinos are Dirac particles and hence there will be no 
neutrinoless double beta decay. If it is indeed found that neutrinos are Majorana particles 
one can think of accommodating that fact into this model by adding some localized new 
sterile neutrinos in the UV brane or some mass term, also localized, for the right handed 
partner of the neutrino. We leave this also for further investigation. 

5 Conclusions 

We have developed a model for ElectroWeak physics embedded in five (non compact) 
dimensions with an exponentially decaying effective metric such that the fifth dimension 
has finite length. The departure of the metric from the AdS case is associated to the 
VEV acquired by some dilaton field <\> for which we have taken (</>) = fi 2 z 2 , z representing 
the fifth coordinate and fi being some mass scale, although the underlying gravity model 
giving room to this particular aspect has not been addressed in the current work. 

In addition, the spacetime has a UV boundary but no IR one. This class of models 
are commonly known as soft wall models in the literature. A smooth fifth coordinate 
dependent Higgs VEV may be justified within these models, turning out in a power law 
behavior. Taking this result into account, we have solved the fermionic equations of motion 
for three different polynomial behaviors of the Higgs VEV: constant, linear and quadratic, 
where the latter is the only one binding the solutions enough for them to be normalizable 
and hence yielding a discrete set of KK excitations. 

We have further solved the corresponding equations of motion for the Higgs and Elec- 
troWeak sectors and found that the smallest eigenmasses are of the same order of mag- 
nitude, both being proportional to fi and thus we are forced to consider the latter at the 
TeV scale. 

Finally, we have found that (the order of magnitude of) the lightest fermionic masses 
behave as a power law of the effective 4d Yukawa couplings, where the exponent is the 
product between the corresponding fermionic 5d bulk mass and the (AdS) radius of cur- 
vature (k). This power law precisely allows us to reproduce the Standard Model hierarchy 
between the fermionic masses with a mild variation for the bulk masses and even order 
1 effective Yukawa couplings. In particular with a bulk 5d mass ranging on ~ (10, 20) \ 
we achieve the correct order of magnitude for the top quark and the neutrinos, respec- 
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tively, with a Planckian curvature, k ~ 10 GeV. While for a much lower curvature, 
k^ 1 ~ 10 4 GeV we reproduce the hierarchy of the masses with 5d bulk masses taking values 
on ~ (2,20) We would like to stress again, however, that our main concern here has 
been to end up with the correct order of magnitude for the SM fermion masses. A more 
realistic study incorporating the CKM matrices and the Electroweak constraints, although 
the main results presented in the current work would not substantially change, alike the 
question concerning the dilaton VEV and its stabilization should be addressed in future 
investigations. 
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A Solution to the Kummer equation for integer first param- 
eter (6) 

The general form of the Kummer differential equation is [10] 

xy" + (b-x)y' -ay = 0. (97) 
For b ^ Z the two linearly independent solutions are 

K=lRMx)s ^±^L^, (98, 

V2 = x 1 - b 1F1 (1 + a- 6,2-6, x) . (99) 
For the confluent hypergeometric function of the first kind, iFi (a, 6, x), we have 

lFl (a, 6, x) « x- b e* + O f i) )+ x- 2 « f - t^- + O ( 1 



,r(a) wy v r (^-«) V^, 

thus a non-exponential asymptotic behavior is found with the linear combination 9 

For 6 G Z the solutions y±, y2 are no longer independent. Actually one easily checks that 

does exist for any n G Z. Indeed, when 6 — > — n with n a positive or zero integer one finds 

W) ( °' 6 ' Z) ~* i^i (1 + a + n, 2 + n, z) , 

and then the function 

U (a, b, z) = ^ ^ (a, 6, z) - ^^^ ^ ^ i*i (1 + a - 6, 2 - M) 
r(l + a-6)r(a) rr , 

= r(i-b)r(6) u ( a M> 

vanishes as b approaches any integer value. Thus two linearly independent solutions are 
given by 

^= S rW lFl(a ' M ' (101) 

V2= J™ { - Z r 1 ^U(a,b,z) = d h U(a,b,z) , (102) 

Sin7T0 b=n 



9 It behaves as a power law. 
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where for the last equality we have used the PHopital rule. In addition notice that 
behaves as a power low when z — > oo since it is the limit of a function proportional to 
U (a, b, z). 

For the particular cases 6 = and 6=1, corresponding to the gauge and the Higgs 
cases, respectively, one easily shows that 

- T (1 + a) zlnz 1F1 (1 + a, 2, z) 



OO 



fc=0 v ; 

— r (a) lnz 1F1 (a, 1, z) , 

with V (x) = {ln[r(x)]} , the di gamma function, and hence we find the (normalizable) 
solutions presented in subsection 3.3. 

B Estimate of the hypergeometric integrals 

To estimate the value of the wave function normalizations and the effective 4d Yukawa 
couplings we will assume To, Co *C 1. Furthermore, since all the integrands are exponen- 
tially suppressed we will approximate the hypergeometric functions to the lowest order 
expansion, disregarding the higher values of the integral variable, that is 



U (a, b, x) 



r(i-6) | r(b-i) xl _ b 



J(l-6 + a) r(a) 



(1 + Or)). (103) 



• H > I 



20 



For the case of the kinetic normalizations we have 
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The leading order integral can be approximated using the L'Hdpital rule when (o — > for 
which we find 
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and then 
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Analogously, for the positive parity spinor ipi we find 
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Notice that Nl (a) = Nr(— a) thus if a were negative the roles of Nl and iV^ would 
simply interchange. 

For the Higgs we take t~ 1 t] (r) ps T ^— lnr and hence 



where we have approximated 
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For the effective Yukawa coupling we start from 
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where we have taken 
■i+a 



6 2 3 



+a U[l--,- + aX)U 



5 1 I 
' 4 ' 2 



<5 2, 



Using again the L'Hopital rule we find 
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and by plugging the normalization constants in the expression (110) and taking zq ~ k, 
the former simplifies to 
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where we have used that Co = ^%\J + v o^ 3 - I n addition, from Eq. (69) we can express 
Co as 

(114) 
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Finally, since we have set /i ~ fov 7 ^ ~TeV we may naturally assume 
and then, from (113) one obtains the corresponding expression in (96). 



• \ < M < 1 

If a were such that 1/2 < |q| < 3/2 we would instead find 
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with ( x ) the digamma function. And finally the effective 4d Yukawa coupling constant 
turn out to be proportional to the physical mass according to 
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C Some technical aspects concerning the solution to the 
fermionic equations of motion 

The system we are concerned with is of the form 



£>i h = g , 
D 2 g = h, 



(116) 
(117) 



where are some (linear) differential operators. The second order (decoupled) system 
is 



D 2 D X h = h, 
D±D 2 g = g. 



(118) 
(119) 



Notice that whenever (116) and (118) hold then (117) and (119) are automatically satisfied 
and vice versa. In addition, let H = { h : R -> R\ D 2 D 1 h = h} and F = { g : R -> R\ D X D 2 g 
be the set of solutions to the second order differential equations (118) and (119), respec- 
tively. Then VlifW (Vj € F), D\h e F (D 2 g G H). These considerations are useful 
when solving the first order constraints. 

Concerning the orthonormality of the solutions to the system (26)-(27) we will show 
that the operators (28)- (29) are indeed hermitian with respect to the usual scalar prod- 



uct 10 , say 
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To see this notice that the second order differential operators can be compactly written as 



OH = H" - 



H 



(120) 



f (M + m D a 3 ) 2 -f'(M + m D a 3 ) - f m' D a 3 

with H = (/i + ,/i_) T , and analogous expression for G = (<?+,<?-) T . Then by means of 
partial integration we find 

(H, OH) = {OH, H) - H^H' + H*H , 

and using the first order constraints 

H' = -(M + m D a 3 )fH-mG, H' = - (M + m D a 3 ) f H - mG , 



(121) 



10 Whenever this scalar product be well defined, which is not the case for the non normalizable solutions. 
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the boundary piece reduces to 

mH ] G -m&H . (122) 

Finally, the boundary conditions (30) allow us to deduce that WG = ffia\{— <J\)G = 
-H^G and thus (H, OH) = (OH, H) for all H, H solving the system (26)-(27). The same 
argument leads the the same conclusion for the operators (29). 

D Fermionic zero-modes profile 

Here we will plot the profiles of the wave functions corresponding to the top quark, for 
which a ~ 13 and 5 ~ 10 -16 




1 2 3 4 5 5 10 15 20 25 30 



Figure 4: Plots for the profiles (from left to right and from top to bottom) h + (blue line), 
g+ (red line), /i_ (green line) and g- (brown line) as a function of C, for a = 13 and 
mk ~ 10~ 17 ; corresponding to the top quark in the Planckian case. 

As can be seen the profiles are peaked towards the UV brane as opposed to the case of 
the hard wall scenario (RS) were the top is peaked in the IR brane. For the case of lighter 
fermions the behavior is qualitatively similar in both scenarios. 
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